Reynolds number effect on the dissipation function in wall-bounded flows 
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The evolution with Reynolds number of the dissipation function, normalized by wall variables, 
is investigated using direct numerical simulation databases for incompressible turbulent Poiseuille 
flow in a plane channel, at friction Reynolds numbers up to Rct ~ 2000. DNS results show that 
the mean part, directly dissipated by the mean flow, reaches a constant value while the turbulent 
part, converted into turbulent kinetic energy before being dissipated, follows a logarithmic law. This 
result shows that the logarithmic law of friction can be obtained without any assumption on the 
mean velocity distribution. The proposed law is in good agreement with experimental results in 
plane-channel and boundary layer flows. 



Despite extensive study, there remain significant ques- 
tions about Reynolds number effects on wall-bounded 
flows. Earlier surveys of data^i^ indicated that Reynolds 
number effects are present in the near-wall region over a 
wide range of Reynolds number. The major difficulty in 
drawing firm conclusions is the accuracy of the measure- 
ments, which invariably are affected by spatial resolution 
and other near-wall measurement issues. 

Direct numerical simulations (DNS) of turbulent flows 
provide detailed turbulence data that are free from such 
experimental ambiguities. During the past two decades, 
the investigations of wall bounded turbulent flows by 
DNS have provided considerable insights into both the 
statistical and structural characteristics of wall bounded 
turbulence. One of the most well-studied turbulent flows 
is the flow in a plane channel, which was simulated by 
Kim et al.^ and by many others since.—"— Moderately 
high Reynolds number simulations have been recently 
performed and the results made available>ii2, 

The aim of the present study is to investigate the 
Reynolds number dependence of the energy dissipation 
function in a turbulent plane channel flow using the re- 
sults of both available DNS databases and our own sim- 
ulations conducted here in order to obtain a wider and 
more complete range of Reynolds number. The simu- 
lation parameters of DNS cases considered are given in 
Tabled 

The present numerical simulations are based on a 
pseudo-spectral code using the Chebychev-tau formula- 
tion in the wall- normal direction (^2) and Fourier ex- 
pansion in the streamwise {xi) and spanwise (xs) direc- 
tions where periodic boundary conditions are appliedi^ 
The number of Fourier/Chebychev modes was selected 
so that the energy spectra are at sufficiently small val- 
ues at large wave numbers, particularly near the wall. 
The flow was driven by a constant streamwise pressure- 
gradient dP/dxi. 

The mean energy dissipation rate per unit volume 0, 
for incompressible flow, is given by>^ 



TABLE 1: Parameters of the turbulent plane channel DNS 
datasets used. Reynolds numbers: R^^ = Ubh/u, Re^ = 

Urh/u. 
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using standard Cartesian tensor notation and summation 
on repeated indices, fi is the dynamic viscosity, Dij and 
dij are respectively the mean and fluctuating part of the 
velocity deformation tensor. The first term on the right- 
hand side represents the part directly dissipated by the 
mean flow ipM, and the second, 4>t, the turbulent part. 
The two terms are usually decomposed as the sum of 
homogeneous and inhomogeneous parts as follows: 
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dxj dxj 



dxi dxj 



and 



dui dui d'^UiUj 

= l^^a ^ — 

oxj oxj oxi oxj 

For fully turbulent flow in a plane channel, the variation 
of mean values in the streamwise and spanwise directions 
are zero. The mean velocity reduces to the streamwise 
component Ui which, like the Reynolds stresses, depends 
only on the wall-normal position X2 ■ Since the gradient of 
the Reynolds stresses is zero at the channel walls, located 
at X2 = ±ft., the dissipation function defined as the 
integral over the channel cross-section of the mean energy 
dissipation rate, is given by 



$ = 



dUi 

dX2 



dui dui 
dxj dxj 



dxo 




FIG. 1: Evolution with the friction Reynolds number Re^ of 
the mean shear "I>^^ (open symbols) and turbulent "3?^ (filled 
symbols) contributions to the nondimensionalized dissipation 
function. , [Eq. Jail]. Symbols in Table II 



For the same reason, the turbulent kinetic energy equa- 
tion (Ref. 11, Eq. 5) integrated across the channel, shows 
that the turbulent contribution to the dissipation func- 
tion is equal to the integral of the turbulent kinetic energy 
production: 



— ^ — da;2 
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where p is the fluid density and uiU2 the Reynolds shear 
stress, and finally 
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Using the integrated streamwise mean momentum equa- 
tion: 



dUi 

dxi 



pUiU2 



X2 



where r„ is the mean wall-shear stress, the dissipation 
function can now be easily evaluated: 



-h 



X2 dUx 
— - — dX2. 
h dx2 



The integration by parts with the no-slip conditions at 
the walls leads to 



Note that an identical relation can be obtained for the 
turbulent flow in a circular pipe (see for example Ref. 
[H ). while for the turbulent boundary layer on a flat plate, 
the dissipation function is related to the streamwise vari- 
ation of the mean and turbulent kinetic energy (see Eq. 
3.15 in Ref. [ij). 

In a dimensionless form, the dissipation function is now 
given by 



dUx 
dxt 



U\U2 



,+ di7i" 



dx\ 



dxl = 2U+, 



where the superscript (-I-) denotes normalization by the 
friction velocity Ur = \J T^j p and the kinematic viscosity 



Figure [T] shows, in a semilogarithmic plot, the evolu- 
tion with i?e^ of and <i>J, the mean and turbuleni 
contributions, respectively. This figure highlights that 
for i?e^ > 500, when $J becomes greater than $ 



reaches a constant value and follows a well defined 
logarithmic evolution. The DNS results give the two re- 
lations 
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5.2 In + 18.27 
512 
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+h 



Uidx2 = 2TwUb 



where Ub is the bulk velocity. This is the classical relation 
for the loss of power in a duct, which is equal to the 
product of the mean pressure gradient —dP/dxi — r^/h 
by the flow rate 2hUb- 



These results are not too surprising according to the 
profiles of the squared mean velocity gradient and the 
mean production for R^^ > 500 shown in figure [21 In 



the region y"*" = 
versal behavior. 



h+ 



^ 30 the profiles exhibit uni- 



Above this position, the square of the 



mean velocity gradient decreases faster than y"*" , giving 
a negligible contribution to the integral, while the mean 
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FIG. 3: Reynolds number evolution of [7^, the bulk veloc- 
ity normalized by friction velocity. Experimental data: (+), 
Zanoun et al. (Ref. HI); (x)^akken et al. (Ref. HI); (■) 

DNS results of Spalart (Ref. E3); , [Eq. (Ell)]. Other 

symbols are in Table [T] 



production decay is close to in a region whose extent 
increases with Reynolds number. 

The skin friction coefficient C/, based on the bulk ve- 
locity is given by: 




2.05 



(0.2) 



This finding shows that for i?e^ > 500, the dissipation 
function is free from low Reynolds number effects and 
corroborates the classical logarithmic law of friction with- 
out the assumption on the mean velocity profile used in 
the classical Prandtl-von Karman analysiS )^^'^'^ based on 
the logarithmic law of the wall 

Ut = -lny++A 

and the defect law 

K h 

where Uc is the mean centerline velocity. The logarith- 
mic skin friction law is obtained by assuming very large 
Reynolds number and combining the two equations (see 
Ref. [H page 573): 



and might be extrapolated to arbitrarily large Reynolds 
numbers. This logarithmic law is compared in Fig. [3] to 
the channel- flow experimental data of Zanoun et al.,^^ 
and Bakken et al}~. Agreement with the data is within 
±1% for i?e, > 500. 



U+ = -lni?e 



A 



1 

K 



Hence, the factor 2.6 in relation (j0.2[) corresponds to a 
von Karman constant k — 0.385 which is close to the ex- 
perimental value of 0.38 obtained by Osterlund et ali^ in 
a zero pressure-gradient turbulent boundary layer. How- 
ever, a small discrepancy with experiment is found for the 
value of the additive constant A since relation (|0.2p leads 
to ^ = 4.65 while Osterlund et al}^ obtained A = 4.1. 

Another noteworthy feature is the good agreement 
of the skin-friction law obtained from Eq. (j0.2p with 
the boundary layer skin-friction coefficient based on 
Ub = (1 — Si/S) Uoo, where Uoo, S and Si are respectively 
the free-stream velocity, the boundary-layer thickness 
and the displacement thickness. This is obvious from 
Fig. U where the experimental data of Osterlund et alt^ 
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and the DNS results of Spalart^ are compared to the 
present relation. This finding requires further analysis 
since the skin-friction coefficient in this flow is given by 
the streamwise mean- momentum variation. 

To summarize, the logarithmic law of friction for plane 
channel flow is the direct consequence of the logarithmic 
evolution of the dissipation function and more precisely of 
the turbulent part. The proposed law of friction is found 
to be in good accordance with boundary layer results. 
A more detailed study including numerical and accurate 
measurements would be useful, particularly for Poiseuille 



pipe flow. 
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